Superfluidity vs. localization in bulk ^He at zero temperature 
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We present a zero-temperature quantum Monte Carlo calculation of liquid *He immersed in an 
array of confining potentials. These external potentials are centered in the lattice sites of a fee 
solid geometry and, by modifying their well depth and range, the system evolves from a liquid 
phase towards a progressively localized system which mimics a solid phase. The superfluid density 
decreases with increasing order, reaching a value Ps/p = 0.079(16) when the Lindemann's ratio of 
the model equals the experimental value for solid *He. 
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I. INTRODUCTION 

The counterintuitive coexistence of superfluidity and 
solid order has been the object of theoretical and exper- 
imental debate for a long time. The recent experimental 
findings by Kim and Chan of non-classical rotational in- 
ertia (NCRI) in solid ''He, both in bulk^ and m a con- 
fined environment,^'^ suppose a real breakthrough in the 
field. Using a torsional oscillator setup and ultrahigh- 
purity ^He, a superfluid density Ps/p ~ 1 — 2% has been 
measured below a temperature T = 230mK. The transi- 
tion temperature from normal solid to supersolid is nearly 
constant with the pressure and no signal of the transition 
is observed in the specific heati^ On the other hand, other 
manifestations of superfluidity like pressure-induced flow 
has not been detected in a recent experiment where solid 
''He was confined in Vycori^ 

On the theoretical side, the pioneering works by An- 
dreev and Lifshitz^ and Chester^ suggested for the first 
time the possibility of Bose-Einstein condensation of va- 
cancies or interstitials in solid ^Hc. Recently, Prokof'ev 
and SvistunoAiS have worked out on the same idea and 
have established that the existence of defects in the solid 
is a necessary condition for observing superfiuidity, i.e., 
the supersolid phase. According to this mechanism for 
superfluidity in the solid phase, ^He could present at very 
low temperature the form of an incommensurate crystal, 
i.e., a quantum solid where the number of lattice sites 
is different from the number of atomsi^ Also, the pos- 
sibility of a disordered solid phase (glass) has been sug- 
gested as an alternative for observing superfluid signals. '^'^ 
However, the hypothesis of incommensurability contra- 
dicts well established experimental data about the abso- 
lutely negligible concentration of vacancies at tempera- 
tures where supersolid has been observedtii 

Following Legett's argumentSfi^ a second plausible sce- 
nario for a supersolid phase would be the one of a com- 
mensurate crystal but with enough exchanges between 
atoms to generate NCRL This fascinating possibility can 
only be thought in a quantum solid such as ^He where the 
displacements of the atoms around the lattice sites are 
incredibly large from a classical point of view. However, 
recent path integral Monte Carlo (PIMC) simulations of 



solid ^He near the experimental transition temperature 
have not detected neither superfluid density^^ nor off- 
diagonal long-range order. ^''^ It is worth mentioning that 
the estimation of both magnitudes is rather difficult using 
PIMC since the temperature is very low and the signal 
is expected to be very small. 

Zero-temperature properties of solid ^He such as its 
equation of state, static structure factor or Lindemann's 
ratio have been calculated using diffusion Monte Carlo 
(DMC)i5 and Green's function Monte Carlo (GFMC).i& 
The results obtained show an excellent agreement with 
experimental data due to both the exactness of the 
method and the accuracy of the interatomic potential 
between helium atoms. However, these calculations use 
as importance sampling trial wave function a Nosanow- 
Jastrow (NJ) model,--^ which is not symmetric under the 
exchange of particles. This is not a drawback for calcula- 
tions of the energy or the spatial structure since the ex- 
change energy has been estimated to be negligible, about 
3/iK/atom,'^'^ but the lack of a correct quantum symme- 
try does not allow for estimations of the one-body density 
matrix and the superfiuid density. Two different models 
have been proposed to overcome this limitation. The 
first one is inspired in the Bloch function of band theory 
currently used for electrons; it is effectively a symmet- 
ric trial wave function but has proven to be worse than 
the NJ model at the variational levelii^ The second one 
corresponds to the shadow wave functionji^ which allows 
for a good variational description of the solid, but not 
yet implemented in a DMC calculation. 

In the present work, we present a microscopic study by 
means of the DMC method of bulk ''He immersed in a lat- 
tice of confining potentials. By increasing progressively 
the strength of these external potentials we can get rel- 
evant information on the interplay between localization 
and superfluidity. Measuring the degree of localization 
through the Lindemann's ratio 7 one can compare the 
physics of this structured phase with the solid one. At 
the experimental value for 7 near the melting point the 
present model shows a non-zero superfiuid fraction, com- 
patible with the existence of a supersolid *He phase at 
zero temperature. 

In the following section, the method used for the calcu- 
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lation is briefly described. Sec. Ill comprises the results 
obtained, and finally, Sec. IV summarizes the work and 
presents the main conclusions. 



II. METHOD 

The iV-body system under study is described by the 
Hamiltonian 

^2 N N N 

^ - E + E n^^^) + E ^cdri - Rj I) , (1) 

1=1 i<j i,j 

with V{r) the interatomic potential between He atoms, 
and 

Vcir) = -Vo exp(-ar2) (2) 

are potential wells centered in the positions Rj, which 
correspond to the sites of a solid lattice. The ground- 
state energy and other relevant properties of the system 
at zero temperature can be calculated exactly, within sta- 
tistical errors, using the diffusion Monte Carlo (DMC) 
methodi^S With DMC, one is able to solve stochasti- 
cally the A'^-body imaginary-time Schrodinger equation 
by performing a random walk of the walkers (R = 
{ri, . . . , tn}), which define the wave function, according 
to a short-time approximation for the Green's function. 

The simulation has been carried out by placing the 
external potentials Vc{r) on the sites of an fee lattice 
and with a number of particles (iV = 108) commensu- 
rate with the lattice. At T = 0, ''He crystallizes in the 
hep lattice but the differences with the fee geometry are 
known to be very small and would not modify the con- 
clusions of the present work. The interatomic potential 
V{r) is the accurate HFD-B(HE) Aziz model2i and the 
trial wave function used for importance sampling is a 
Jastrow factor V' = exp[^^^- w(r.y )], with a McMillan 
correlation u{r) = — 0.5(6/r)^. The parameter b = 1.3cr 
(<7 = 2.556A) is the same for all the simulations and the 
attractive potentials on the sites J^jl are varied on its 
strength and range among the values Vq = (500 — 800)K 
and a = (1 — 2)cr^^. These very sharp potentials are 
required to confine the atoms around the sites at a quan- 
titative level comparable with the one of solid ''He. It 
is worth noticing that the localization arises due only 
to the branching term in the DMC algorithm since we 
have chosen not to include confining terms in the trial 
wave function to ensure a fully symmetric wave function. 
In order to reduce any residual bias in the calculation a 
preliminary study on the time-step and on the number 
of walkers has been carried out. In the range of confine- 
ment simulated, the energies are well behaved and their 
statistical errors are smaller than 1%. 



III. RESULTS 

The influence of confinement on the system has been 
studied at a density p = 0.491cr~'^, which corresponds 
to an experimental pressure of ~ 34 bar, close to the 
melting point of solid "^He. In Fig. 1, we show a projec- 
tion on the X — y plane of imaginary-time trajectories of 
some walkers and for a given evolution time for both the 
present model and a DMC simulation using the NJ trial 
wave function (DMC-NJ) at the same localization degree. 
As one can see, the displacements of the atoms around 
the lattice sites are in both cases rather large according 
to the quantum nature of "^He. The difference between 
both simulations relies on the exchanges between atoms 
at neighboring sites, that are observed in the confining 
model, and which are forbidden by construction in the 
NJ description. With this configuration of the confining 
potentials the probability of hopping, i.e., the quotient 
between the number of MC steps where a hop is pro- 
duced and the total number of steps along a simulation, 
is - 3%. 

In order to characterize the degree of localization as 
a function of the range and strength of Vc{r) we have 
calculated the density profile /i(r) of ''He atoms around 
the lattice sites {Rj}. From them, one can estimate the 
mean squared displacement 

/>oo 

{u^} =An niry^dr , (3) 
Jo 

and then the Lindemann's ratio 7 = ^/{u?)/a, a being 
the nearest-neighbor distance in the lattice. In Fig. 2, 
results for /x(r) for three confining models are plotted in 
comparison with the DMC-NJ result of sohd "He. The 
values for 7 of the three cases are 0.192(1), 0.264(2), and 
0.294(3), to be compared with the DMC-NJ value 7 = 
0.257 which agrees with the experimental value 7'^'^P* = 
0.26.22 AH the fines in Fig. 2 are fits to the DMC results 
using a sum of two Gaussians, a simple model which has 
proved to reproduce accurately the computed data. The 
figure shows that the present model is flexible enough 
to simulate from a strongly confined system, with 7 = 
0.192(1), until a more delocalized one, with 7 = 0.294(3), 
crossing the value of the Lindemann's ratio of solid '*He. 
As one can see, the curve of the model for 7 — 0.264(2) 
in nearly indistinguishable from the DMC-NJ result. 

The main objective of the present work is the estima- 
tion of the superfluid density as a function of the de- 
gree of localization of particles around the lattice sites. 
The superfluid density of a bosonic system may be com- 
puted with DMC by extending the winding-number tech- 
nique, originally developed for PIMC calculations, to zero 
temperaturei^ Essentially, this is obtained by calculat- 
ing the diffusion constant of the center of mass of the 
particles in the simulation box for an infinite imaginary 
time 

^ = lim ^ (^) , (4) 
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FIG. 2: Density profile, /i(r), of tlie ''He atoms around the site 
positions {Rj}. Tlie solid line corresponds to the DMC-NJ 
model for solid *He (7 = 0.257). Dotted, dashed, and long- 
dashed lines stand for the present model using potentials Vc(r) 
that yield 7 = 0.192(1), 0.264(2) and 0.294(3), respectively. 



where Ds(t) = ((RcM(r)-RcM(0))2), andDo = rc l^m. 
Results obtained for ps corresponding to several config- 
urations of the lattice potentials K are plotted in Fig. 
3. In this figure, Psj p is shown as a function of the Lin- 
demann's ratio of the same configuration and for values 
0.18 < 7 < 0.30. The results for /0s(7) are well parame- 
terized by the function 



1 + cxp 



7-70 



(5) 



shown as a line in Fig. 3. The parameters in Eq. © 
are 70 = 0.313(4) and uj = 0.021(1). From Eq. © one 
can estimate the prediction of our model for the range of 
more interesting values for 7. In particular, considering 
^expt ^ Q the superfluid fraction is p^/ p = 0.079(16). 
It is worth noticing that the present results for Ps/p(7) 
are qualitatively similar to the ones obtained by Saslow^"* 



FIG. 3: The superfluid fraction as a function of the Linde- 
mann's ratio 7 in the localized model. The line is a fit @ to 
the computed data. 



assuming an uncorrelated model where all particles de- 
velop a common phase function. However, his most re- 
cent calculations predict smaller values for the superfluid 
fraction than ours, and in particular, for solid '^He a value 
pjp = 0.022. 

The connection between the localized model here pre- 
sented and the NJ description of solid ''He has been made 
so far by comparing the Lindemann's ratio of both calcu- 
lations. To make this comparison more complete, we have 
explored if there are other structure properties which can 
also be comparable for similar values of 7. To this end, we 
have computed the two-body radial distribution function 
g(r), which is proportional to the probability of finding 
two particles a distance r apart, and its Fourier trans- 
form, the static structure factor S{k). 

In Fig. 4, results of g{r) for two different configura- 
tions of the confining potentials Vc{r) are shown in com- 
parison with the DMC-NJ result at the same density. 
These functions, as also the density profiles p{r) shown 
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FIG. 4: Two-body radial distribution function of the localized 
model corresponding to superfluid fractions ps/ p = 0.08 (long- 
dashed line) and 0.035 ( short-dashed line). The solid line 
stands for the DMC-NJ result at the same density. 




FIG. 5; Static structure factor of the localized model cor- 
responding to Pa/p = 0.11 and 7 = 0.264 (squares). The 
DMC-NJ result is also shown for comparison (triangles). 



in Fig. 2, have been calculated using a pure estimation 
technique to avoid the residual bias coming from the trial 
wave function that the direct simulation output (mixed 
estimation) presents. As one can observe in the figure, 
when the Lindemann's ratio of the confining model is 
very close to the one of the NJ model both curves co- 
incide. In the same figure, a result corresponding to a 
smaller 7 value shows clearly a more pronounced peak 
reflecting its higher localization around the lattice sites. 

The solid-like character of a given configuration is 
much more clear by looking at the static structure func- 
tion S{k) = l/Ar(p_kPk), with pk - E^ie'""""'- than 
g{r). In Fig. 5, results of S{k) for the present model are 
compared with the DMC-NJ result for solid '*He. The 
localized model result shown in the figure corresponds to 
a Lindemann's ratio 7 = 0.264(2), close to the DMC- 



NJ value for 7, and a superfiuid fraction ps/p = 0.11(1). 
Both results show the presence of the two main first peaks 
at the same k values, characteristic of the underlying fee 
lattice. The strength of the peaks of the localized model 
is slightly lower than the one of the DMC-NJ calculation 
since the values for 7 in both models are not exactly the 
same, but it is clear that the present model is able to 
reproduce the expected trend of a solid phase in S{k). 

A slight trend of increasing Ps/ p with pressure is ev- 
idenced on Kim and Chan's experiments,^ contrarily to 
what it occurs in superfluid ^He; the authors conjecture 
that this may be due to the presence of grain boundaries 
formed in the growth of the crystal and which affect the 
coherence in the superflow. Regarding this issue, we have 
computed Ps/p at higher density, pu = 0.529ct~^, and 
external constraint yielding 7 = 0.264(2) and Ps/p ~ 
0.13(3) at density p = 0.491(7"^. The pressure of bulk 
solid ''He at density p^ (~ 64 bars) is still on the range 
of supersolid experiments. We flnd ps/p = 0.10(2) and 
7 = 0.239(5) (7 = 0.246(3) for the DMC-NJ simulation 
at Pu). Therefore, within the statistical uncertainties of 
the model we do not observe pressure dependence oi Ps/p 
in the range of pressures analyzed. 



IV. CONCLUSIONS 

We have studied the competing effects of localization 
and superfluidity of bulk ^He immersed in a lattice of 
external conflning potentials. The simulations have been 
carried out at zero temperature using the DMC method. 
From the curve ps/p(7), and using the experimental value 
of 7 of solid ^He we obtain an estimation of the superfluid 
fraction ps/p — 0.079(16), a value significantly larger 
than the experimental measure {pf^^^ / p = 0.01 — 0.02). 
On top of the possible drawbacks of the localized model 
for describing the ^He solid phase, there are two effects 
that could help to explain the difference between the ex- 
perimental measure and the present prediction. They 
are, first, the temperature (the simulation is at zero tem- 
perature), and second, the possible infiuence of impuri- 
ties in the experiments. The latter can be specially rele- 
vant since, as commented by Kim and Chan,^'^ a crucial 
point to observe superfluidity is the as much as possible 
elimination of impurities in the sample. The influence on 
the superfluid signal of '^He impurities and/or incommen- 
surability can be also analyzed in the present framework. 
Work is in progress along this direction. 
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